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PREFACE 


This  thesis  is  a  continuation  of  the  work  done  by  Yeakel 
in  the  control  of  linear  periodic  systems.  While  specific 
cases  are  looked  at  in  this  study,  the  techniques  are 
applicable  to  linear  periodic  systems  in  general.  This 
field  is  Just  beginning  to  be  explored  in  depth.  With  the 
advent  of  computers  to  perform  many  of  the  tedious  numerical 
integrations  the  problem  becomes  workable.  With  this  thesis 
I  hope  to  add  one  more  building  block  to  this  field.  More 
work  needs  to  be  done  before  we  can  fully  understand  the 
nature  of  periodic  systems,  and  how  to  control  them.  By 
filling  in  all  of  the  bits  and  pieces,  eventually  the  field 
will  be  conquered. 

Much  of  this  thesis  work  was  done  on  the  computer,  a  tool 
I  was  not  completely  comfortable  with.  I  would  like  to  thank 
all  the  people  who  helped  me  understand  the  operation  of  the 
computer  system,  especially  Dr.  V/illiam  Wiesel,  and  Captain 
Hugh  Briggs.  I  would  also  like  to  thank  the  members  of  my 
faculty  board.  Captain  DeWispelare,  Major  Wallace,  and  Dr. 
Wiesel.  Most  importantly  I  would  like  to  thank  my  adviser 
Dr.  Robert  Calico  for  the  guidance  and  insight  he  has  given 
me  this  past  year. 
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ABSTRACT 


'-^The  attitude  of  a  spinning  symmetric  satellite  in  an 
elliptical  orbit  was  analyzed.  The  linearized  equations 
were  formed  from  which  the  stability  was  determined  using 
Floquet  theory.  The  majority  of  satellite  configurations 
examined  exhibited  two  unstable  modes.  Control  laws  using 
pole  placement  techniques  were  implemented  which  stabilized 
the  satellite.  A  scaler  control  provided  stability  but  did 
not  allow  for  exact  pole  placement  of  the  controlled  system. 
Multiple  input  control  gave  exact  pole  placement  but  required 
complicated  time  dependent  gains. 
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CHAPTER  1 


INTRODUCTION  AND  BACKGROUND 

Since  the  beginning  of  space  flight  in  the  late  1950 's 
designers  have  been  concerned  with  maintaining  the  attitude 
of  their  satellite  with  respect  to  some  fixed  reference. 

Today  this  is  still  a  major  concern,  as  the  use  of  earth 
orbiting  satellites  is  vital  not  only  for  weather  and 
communications  problems,  but  also  for  national  security 
interests.  The  need  to  keep  an  unmanned  satellite  "looking" 
at  a  particular  point  on  the  earth’s  surface,  whether  it  be 
a  storm  front  or  a  foreign  military  installation,  is  an 
important  problem  facing  satellite  designers.  Various  methods 
have  been  used  in  the  past  to  achieve  this  "line  of  sight" 
precision,  involving  both  passive  and  active  control  mech¬ 
anisms.  For  satellites  in  nearly  circular  orbits  a  gravity 
gradient  stabilization  has  been  used.  This  method  has  many 
restrictions  on  satellite  design  and  also  requires  nearly 
circular  orbits.  Gravity  gradient  stabilization  yields  very 
low  natural  frequencies,  and  hence  any  disturbances  take  a 
long  time  to  die  out;  not  a  good  condition  for  precision 
serveillence .  More  active  attitude  control  devises  are  used 
such  as  magnetic  torquers,  mass  movement  and  momentum  exchange 
devises,  and  Jet  thrusters.  Each  has  its  own  advantages  and 
disadvantages  depending  on  the  mission  reqirements.  The  most 
important  factor  in  selecting  a  control  devise,  and  for  the 
entire  satellite  package,  is  weight.  The  lighter  the  overall 
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system  is,  the  higher  orbit  it  can  be  put  into,  and  the  lonper 
it  will  stay  in  space.  Therefore  the  lighter  the  control 
system  the  more  weight  is  available  for  payload.  The  more 
payload  implies  a  higher  degree  of  mission  accomplishment  or 
sophistication . 

One  active  control  devise  long  employed  is  the  gyroscope. 

A  gyroscope,  having  a  large  amount  of  angular  momentum  with 
torque  free  motion,  results  in  a  nearly  constant  angular 
momentum  vector.  Thus  it  is  ideal  for  an  attitude  reference 
devise.  Gyros  have  long  been  used  for  aircraft  navigation 
and  so  their  use  in  spacecraft  seemed  natural.  To  simplify 
the  problem  the  idea  became  to  make  the  entire  satellite 
gyroscopic  by  spinning  it  about  some  axis.  As  early  as  1963 
Kane  and  Shippy  investigated  the  stability  of  a  spinning 
satellite  in  a  circular  orbit,  (5;111).  This  yielded  equations 
which  were  nonlinear  and  nonautonomous .  Later  in  1966  Kane 
and  Barba  investigated  the  problem  of  a  symmetric  satellite 
in  an  elliptical  orbit,  spinning  about  its  axis  of  inertial 
symmetry,  (4;402).  Again  the  system  was  nonlinear  and  non¬ 
autonomous  . 

In  both  cases  the  equations  were  linearized  about  an 
equilibrium  point  giving  linear  equations  with  periodic 
coefficients  whose  stability  was  checked  using  Floquet 
theory.  In  the  case  of  a  satellite  in  an  elliptical  orbit 
the  stability  was  found  to  be  dependent  on  the  inertia  prop¬ 
erties,  the  orbit  eccentricity,  and  the  spin  rate  of  the 


satellite , 


As  the  need  for  more  precise  line  of  si&ht  capability  has 


increased,  the  use  of  gravity  gradient  stabilization  and 
other  schemes  overly  restrictive  of  satellite  design,  has 
become  undesireable i  Some  sort  of  on  board  control  system, 
or  autopilot  is  needed.  For  the  case  of  the  symmetric  spin¬ 
ning  satellite  in  elliptical  orbit,  many  of  the  satellite 
configurations  are  inherently  unstable.  As  the  eccentricity 
is  increased  more  of  the  cases  become  unstable.  A  control 
scheme  needs  to  be  implemented  so  that  these  unstable  con¬ 
figurations  ean  be  made  stable  and  thus  operational.  The 
linearized  system  of  equations  describing  this  case  (Kane  and 
Barba)  have  periodic  coefficients  and  so  classical  control 
techniques  used  for  systems  with  constant  coefficients  are 
not  directly  applicable  and  must  be  modified.  Very  little 
work  has  been  done  on  control  of  periodic  systems.  Shelton 
developed  a  control  scheme  for  periodic  systems,  in  his  thesis, 
for  a  satellite  in  orbit  about  the  earth  moon  Lagrange  point 
L4.  (Ref; 9)  Yeakel  used  the  same  principles  in  developing  a 
control  scheme  for  the  case  of  an  unsymmetric  satellite  in  a 
circular  orbit,  and  a  symmetric  satellite  in  an  elliptical 
orbit.  (Ref;ll) 

This  control  scheme  involves  transforming  the  state  variables 
to  modal  variables  so  the  designer  can  look  at  the  individual 
modes  and  their  respective  stability.  Control  terms  can  then 
be  applied  to  the  individual  modes.  The  control  terms  are 
based  on  feedback  of  the  selected  modes,  modified  by  selected 
gains.  In  Shelton's  and  Yeakel's  work  the  controller  v;as 
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designed  to  stabilize  a  system  where  one  mode  was  unstable. 
Many  of  the  oases  presented  in  this  study  have  two  unstable 
modes.  This  makes  the  control  a  more  ■-implicated  process  due 
to  the  coupling  of  the  unstable  modes.  This  results  in 
unpredictable  results  using  the  scaler  control  developed  by 
Shelton  and  Yeakel.  The  use  of  multiple  input  control  for 
a  system  with  two  unstable  modes  is  therefore  developed.  This 
control  scheme  gives  predictable  results,  but  has  its  own 
special  drawbacks,  namely  complicated  gains. 

This  study  addresses  the  case  of  a  symmetric  satellite  in 
an  elliptical  orbit,  spinning  about  its  axis  of  inertial 
symmetry,  where  two  modes  are  unstable.  The  stability  is 
dependent  on  the  eccentricity,  inertia  properties,  and  the 
spin  rate  of  the  satellite,  and  the  majority  of  cases  involve 
two  unstable  modes.  Thus  the  ability  to  adequately  control 
such  a  case  is  extremely  important.  It  is  the  purpose  of 
this  study  to  investigate  both  scaler  control  schemes  and 
multiple  input  control  schemes  in  stabilizing  such  a  system 
with  two  unstable  modes. 


CHAPTER  2 
THEORY 


In  this  chapter  the  equations  of  motion  of  an  inertially 
symmetric  satellite  in  an  elliptical  orbit  about  the  earth 
are  developed.  The  satellite  is  spinning  about  its  axis  of 
inertial  symmetry.  Also  the  equations  describing  the  attitude 
motion  of  the  satellite  are  presented.  These  attitude  equa¬ 
tions  will  be  linearized  about  an  equilibrium  point  to  des¬ 
cribe  the  system  in  a  conventional  linear  form.  The  system 
is  time  periodic  due  to  the  nature  of  the  orbital  motion, 
therefore  Floquet  theory  can  be  used  to  make  a  statement 
about  the  stability  of  the  attitude  motion.  An  unstable 
system  can  be  made  stable  by  inserting  a  feedback  control  term 
in  the  model.  This  control  term  can  either  be  scaler  control 
or  multiple  input  control.  Both  cases  are  developed. 

EQUATIONS  OF  MOTION 


The  development  that  follows  is  taken  from  Kane  and  Barba 
(4;402).  Two  basic  equations  describe  the  trajectory  of  a 
satellite  in  orbit  about  a  spherically  symmetric  attracting 
body.  These  equations  are  developed  in  any  introductory 
astrodynamics  text  (10;31).  The  equations  are; 


-  2  3 

..  .2  n  a  « 

r  -  ri*  +  —  =  0 

r 


(2.1) 
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and, 


r^i  =  a^n  >/ 1  -  (2.2) 

where 

n  =  2t/T  (2.3) 

and  where  T  is  the  period  of  the  orbit.  The  orbital  elements 
r,i/,a,  and  e  are  shown  in  Figure  2.1; 


Figure  2.1  Elliptical  Orbit  Elements 
The  variables  may  be  nondlmensionalized  by  changing  to  a 
new  system  such  that: 


T  =  nt 

By  rearranging  equation  (2.2), 


(2.4) 


V  = 


i^n  >/l  -  e^ 


substituting  this  into  equation  (2.1)  gives; 


(2.5) 


4„2,,  2.  23 

r  -  - - 5 -  +  — ^  -  0 

r  r 


(2.6) 
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Now  using  the  nondimensional  variables  and  denoting  diff¬ 
erentiation  with  respect  to  t  by  primes: 


2 

ft  I  +  ig-  ^  +  ^2  “  ° 

The  initial  conditions  for  this  differential  equation  are 
found  by  assuming  the  satellite  starts  at  perigee,  therefore; 


or, 


rp  =  a(l  -  e)  Tp  =  0 

f(0)  =  1  -  e  r'(0)  =  0 


(2.8) 


Nondimensionalizing  equation  (2.5)  gives; 


and, 


•  t  -2  y/l  -  e^  V 


(2.9) 


(2.10) 


It  is  worth  noting  that  f  is  a  periodic  function  of  r . 

Nov/  the  equations  to  describe  the  attitude  motion  of  the 
satellite  can  be  developed.  First  an  orbital  reference 
frame  A  is  defined.  Frame  A  is  such  that  points  outward 
along  the  orbit  radius  r,  Ag  is  perpendicular  to  A^^  in  the 
orbital  plane  such  that  at  perigee  the  satellite  is  moving 
only  in  tne  positive  Ag  direction,  and  Ag  is  perpendicular 
to  both  Aj^  and  Ag,  pointing  out  of  the  orbital  plane.  A  body 
axis  frame  X  is  obtained  by  meams  of  a  1-2-3  rotation  through 
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the  angles  9^,  6^,  and  9^ 


Since  the  satellite  is  assumed  to 


be  inertially  symmetric  a  nodal  axis  system  C  is  introduced 
which  is  obtained  from  A  by  means  of  a  1-2  rotation  through 
9^  and  tfg.  See  Figure  2.2. 


Figure  2.2  1-2-3  Rotation 

Expressing  the  inertial  angular  velocity  components  in  the  C 
frame,  including  the  orbital  rate  •'A^; 

«1  *  fl^costf^  ~  I' cosffj^sintf^ 

«2  =  ^2  •'Sin^i  (2.11) 

=  9^sin9^  +  9^  +  ^cos9^cos9^ 

The  moment  equation  (7;468)  gives; 

with  the  applied  moments  being  (4;403); 


0 


Ml  = 


=  3^(1 

M3  =  0 


J)  3111^2003^2 


(2.13) 


Now  from  Euler's  equations; 


( 

)w2W3 

M. 

“1 

+ 

T 

“2 

+ 

( 

T 

_j 

^“l"3 

“3 

+ 

( 

T 

)wiw  2 

(2.14) 


where  the  a's  are  the  angular  accelerations  and  the  I's  and 
J's  are  defined  from  the  mass  moment  of  inertia  matrix; 


0 

0 

J 


(2.15) 


for  a  symmetric  satellite.  Defining  an  inertia  parameter  K 
such  that; 

K  =  -  (2.16) 


The  inertia  properties  can  be  fully  described  in  this  one 
parameter.  The  acceleration  terms  Oi  are  found  using  the 
equation; 


a  = 

dt 


X  (2.17) 


which  given  in  component  form  are; 
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=  -VcosdjSintf2+  V  -  tfgcos  S^^cos +  9^  sinff^) 

+  e'j^cose^  -  e^e^sine^  + 


ttg  =  i''sin®^  +  +  «j^cos  «^sin«2)  +<>2  “  ®3  *1 


(2.18) 


ycos9^cos9^  -  i’tfj^sintf^cos  #2  ~  »<«2Costf^sin  ^2  + 
9^  sln9^  +  tfj^^2COS»2  +  ^3 


Now  having  solved  for  the  a's,  equation  (2.14)  becomes; 


+  Ka»^a»2  ~  ® 

2  3 

«2  “  ^"i^a  “  Ksinff2cos02/  f 


(2.19) 


«3 


=  0 


Consider  the  motion  such  that  0^^  =  ©2  ®  0  and 
where  =  n(  1  +  a )  with  a  a  constant  representing  the  spin 
rate  about  the  C-  axis.  Under  these  conditions  the  accelera- 

O 

tions  a^=  0  and  the  external  moments  M^=  0.  This  represents 
an  equilibrium  point  with  the  axis  remaining  perpendicular 
to  the  orbital  plane.  Now  if  the  satellite  is  perturbed  from 
this  equilibrium  state,  there  is  no  guarantee  the  motion  will 
return  the  satellite  to  this  orientation.  To  determine  the 
stability  of  the  satellite,  state  variables  are  defined  such 
that ; 


Xi  =  X2  =  02  X3 


x^  =  0^ 


(2.20) 


Linearizing  about  the  equilibrium  point  described  above,  using 
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equations  2.11,  2.18,  and  2.19,  gives  the  perturbation 
equations  in  state  vector  form; 


where 


-  T 

X 


x'  =  (  A j  X 

=  [©11  ©2.  ©{.  ©2) 


(2.21) 


hence 


and 


0 

0 

1 

0 

®1 

►  = 

0 

0 

0 

1 

< 

^2 

0'- 

^31 

0 

^34 

©1 

®2 

V  ^ 

.  t* 

1/ 

■ 

^42 

-A34 

0 

®2- 

(2.22) 


^11  ~  ^12  ~  ^14  “  ^21  *  ^22  '  ^23  “  ^33  “  ^44  ®  ° 
^13  ^  ■^24  ”  ^ 


A31  =  -i/'(e^  +  Ki.'  +  Ke^) 
"^3  2  "^41 


(2,23) 


'34 


=  -  e'  -  K(  I-'  +  0')  =  -A 


43 


^42  “  ^31  “  ^ 


This  expresses  the  A  matrix  in  terms  of  the  variables 
^'and  f  .  The  variables  v'  and  "  can  be  eliminated  from  the 
matrix  by  noting  equations  (2.9)  and  (2.10).  can  be 

eliminated  by  noting  the  third  of  equation  (2.14);  =  0 

which  when  linearized  is; 


jJ  +  0^  =  constant 


(2.24) 
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let  that  constant  be  (  <*  +  1  ) ,  then; 

e '  =  (  1  +  «  )  -  - - ^  (2.25) 

3  f2 

Mal<ing  these  substitutions  in  the  A  matrix,  the  elements  are 
only  functions  of  the  variables  K  and  its  derivative  f. 
Additionally  the  elements  depend  upon  K,  a,  and  e,  specifically; 

■^11  *  ‘'^12  “^14  *  ^21  "^22  ^  ^22  ^  ^22  '^44  ° 

^13  "  ^24  ^ 

A31  =  (  1  -  -  (  1  +  a)(  1  +  K  )(  1  -  e^)'''f^ 

A32  =  -2(  1  -  e2)r”2r'  =  (2.26) 

A34  =  2(  1  -  e^)'^r^  -  (  1  +  a  )(  1  +  K  )  =  -A^^ 

A42  =  -  3Kr=' 


FLOOUET  THEORY 

It  is  important  to  note  the  A  matrix  is  periodic 
with  period  T,  that  is; 


A(t)  =  A(  t  +  T  )  (2.27) 

where  T  is  the  period  of  one  orbit.  The  stability  of  the 
motion  can  now  be  analyzed  using  Floquet  theory.  Floquet 
theory  is  used  because  the  matrix  A  is  a  time  periodic 
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function  and  does  not  lend  itself  to  the  same  stability 
analysis  as  that  of  a  constant  A  matrix.  The  fundamental 
matrix  for  this  system  must  satisfy  the  matrix  differential 
equation ; 


<t'(r)  =  A(t)<Kt)  (2.28) 

This  ♦  is  one  fundamental  matrix  of  the  system.  Floquet  theory 
states  that  the  4>(r)  may  be  written  in  the  form; 

♦  (t)  =  P(t)  (2.29) 

where  P(t)  is  a  periodic  matrix  with  the  same  period  as  the 
original  system,  and  the  matrix  F  is  a  constant  matrix.  It 
is  desireable  to  have  the  F  matrix  in  diagonal  form,  (i.e. 
Jordan  form).  This  may  be  done  by  means  of  a  similarity 
transformation  and  a  new  fundamental  matrix  ^  is  therefore 
given  by; 


^(r)  =  F(r)  e'^’’  (2.30) 

F(t)  is  periodic  such  that, 

F(r)  =  F(  r  +  T)  (2.31) 

and  J  is  the  Jordan  form  of  F  . 

Now  returning  to  the  form  of  equation  (2.29)  to  write; 

*(  r  +  T)=P(r  +  T)e^^’’  +  (2.32) 
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but 


P(r)  =  P(  r  +  T  )  (2.33) 

therefore , 

<*>  (  r  +  T  )  =  P(r)  =  <|.  ( r )  (2.34) 

or 

*(  r  +  T  )  =  *(r)  C  (2.35) 

Here  C  is  a  constant  matrix  called  the  monodromy  matrix.  The 
eigenvalues  of  the  F  matrix  are  called  the  characteristic 
exponents,  .  The  eigenvalues  of  C  are  called  the  character¬ 
istic,  multipliers  .  By  letting  t  =  0  tlie  monodromy  matrix 
can  readily  be  found; 

<t>"^(0)<i>(T)  =  C  (2.36) 

If  is  the  principle  fundamental  matrix; 

♦(0)  =  I  (2.37) 

<f>  can  be  solved  numerically  by  recalling  equation  (2.28) 
and  given  the  initial  conditions  (2.37). 

The  eigenvalues  of  the  monodromy  matrix  ai'e  related  to  the 
characteristic  exponents  by; 


The  stability  of  the  system  is  determined  by  the  character¬ 
istic  exponents.  If  all  of  the  characteristic  exponents  are 
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less  than  zero  the  system  is  stable.  If  Just  one  of  the 
characteristic  exponents  is  greater  than  zero  the  system  is 
unstable.  For  the  case  where  the  characteristic  multipliers 
are  complex,  the  relationship  between  them  and  the  character¬ 
istic  exponents  is; 


^i 


In  r^  +  J  arg  p 


(2.39) 


where 


r . 
1 


2  ^  ,  .2 
+  Im(p^) 


(2.40) 


arg  =  arctan(  Im( p^ ) /Re ( p^ ) ) 

In  this  case  all  real  parts  of  the  characteristic  exponents 
must  be  nonpositive  for  stability.  A  positive  real  part 
means  the  system  is  unstable. 

The  system  stability  can  be  determined  from  the  monodrom.y 
matrix,  but  the  periodic  solution  matrix  will  be  needed  for 
implementation  of  control.  A  solution  to  the  original 
system  can  be  written  as  (3;234); 


^(J-)  =  \('(t)^“^(0)x(0)  (2.41) 

v/here  ^  is  any  fundamental  matrix  and  x(0)  is  the  initial 
condition  vector,  assumed  given.  Using  the  fundamental 
matrix  of  equation  (2.30),  and  taking  its  inverse  gives; 

^-^(r)  =  F-^(r)  (2.42) 


where  at  t  =  0, 


15 


(2.43) 


lA"^(0)  =  e°  F  ^(0)  =  F"^(0) 
using  this  in  equation  (2.41); 

x(t)  =  F(t)  F“^(0)x(0)  (2.44) 

It  is  important  to  note  that  F(  r  +  T  )  =  F(t).  in  terms  of 
the  principle  fundamental  matrix  x(r)  may  be  expressed  as; 

x(r)  =  <t.(r)4,-^(0)x(0)  (2.45) 

Evaluating  at  r  =  T , 

x(T)  =  <l>(T)«i.“^(0)x(0)  (2.46) 

where  from  equation  (2.36), 

<I)(T)  *"^(0)  =  <I»(T)  I  =«I>(T)  =C  (2.47) 

where  C  is  the  monodromy  matrix.  Evaluating  equation  (2.44) 
at  T  =  T  and  then  equating  to  equation  (2.46)  gives; 

F(T)  e'^'^  F“^(0)x(0)  =  Cx(0)  (2.48) 

or  by  rearranging, 

F(T)  e'^'^  -  C  F(0)  =  0  (2.49) 

Since  F(r)  is  periodic,  F(0)  =  F(T)  so  equation  (2.49) 
becomes ; 

F(0)  -  C  F(0)  =  0  (2.50) 
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JT 

The  matrix  e  is  diagonal  and  hence  the  individual  columns 
of  the  F(0)  matrix  can  be  written  as; 

-  C  f^(0)  =  0  (2.51) 

JT 

with  being  the  1 ' th  diagonal  element  of  e  .  Equation 

(2.51)  may  be  recognized  as  a  basic  eigenvalue  problem. 

Written  in  a  more  familiar  form; 

|.I  -  c)  7(0)  .  0  (3.52) 

Hence  the  F(0)  matrix  is  a  matrix  whose  columns  are  the 
eigenvectors  of  the  monodromy  matrix  C. 

Referring  to  equation  (2.30)  again,  taking  the  derivative 
of  ^  with  respect  to  r  gives; 

^*(t)  =  F'(7)  e"’’  +  F(r)  j  e'^’'  (2.53) 

Since  ^  is  a  fundamental  matrix  it  raust  satisfy  equation 
(2.29)  ,  i.e. 

^'(’■)  =  A(t)  ^(r)  (2.54) 

Using  the  expression  for  ^  and  in  equation  (2.54)  gives, 

F'(>-)  e'^’’  +  F(7)  j  e'^  ^  =  A(r)  F(r)  ^  (2.55) 

By  cancelling  out  the  e'^ terms  and  rearranging; 


17 


F'  =  AF  -  FJ  (2.56) 

This  equation  alon?  v;ith  the  initial  condition  matrix  found 
in  equation  (2.52)  can  be  solved  numerically  for  ?(»•). 

Summarizing  the  steps  to  determine  the  stability  of  the 
system; 


X '  =  A(  ’^ )  X 

are  as  follows; 

1)  Determine  the  principle  fundamental  matrix  by  numerically 
integrating  equation  (2.28)  v/ith  initial  conditions 


2) 


<I>(0)  =  I 
4>'  (  »■  )  =  A(  t)«|>(  r  ) 

Determine  the  monodromy  matrix  C  from  equation 


(2.37) 

(2.28) 

(2.36) 


<i>”^(0)<f»(T)  =  c 


(2.36) 


3)  Find  the  eigenvalues  of  the  monodromy  matrix  C, ,  and 
therefore  the  characteristic  exponents  by; 


X 


i 


(2.38) 


If  all  the  X^  have  negative  real  parts  the  system  is  stable, 
otherv/ise  it  is  unstable. 

Summarizinor  the  steps  to  determine  the  periodic  F  matrix 
of  equation  (2.30)  are  as  follows; 
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1)  Find  eigenvectors  of  the  monodromy  matrix  from  equation 
(2.52),  These  form  the  initial  conditions  for  the  F  matrix. 


[p  I  -  C  ]f(0)  =  0  (2.52) 

2)  Solve  equation  (2,56)  for  F  by  numerically  integrating. 

F'  =  AF  -  FJ  (2.56) 

Now  all  elements  of  the  general  solution  form  of  equation 
(2.30)  are  known. 


i{r)  =  F(  O  '  (2.3C) 

where  J  is  the  Jordan  matrix,  (7; 267). 

Equations  (2.28)  and  (2.56)  can  be  solved  num.erically , 
given  the  initial  conditions.  The  elements  of  the  F  matrix 
are  needed  as  functions  of  r  and  since  they  are  periodic,  may 
be  stored  as  Fourier  coefficients. 


CONTROL  THEORY 


The  modal  control  of  linear  constant  coefficient  systems 
is  well  understood.  The  application  of  modal  control  for 
linear  periodic  systems  is  however  not  well  developed  and  is 
the  subject  of  this  section.  Consider  a  change  of  state 
variables  from  x  to  a  new  modal  variable  v,  given  by; 


x(  r)  =  F(  r  )■?(  r  ) 


(2.57) 
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where  F(t)  is  the  same  matrix  as  in  equation  (2.30).  Taking 
the  derivative  of  x  with  respect  to  t  gives; 


X'  =  F'i>  +  Fn'  (2.58) 

Substituting  equation  (2.57)  and  (2.58)  into  the  original 
system  of  equations  (2.21); 


F'  »»  +  F>» '  =  AFn 


By  rearranging, 

1  »  =  F'^AF'^  -  F'^F'ir 

Using  equation  (2.56)  to  change  equation  (2.60)  to; 


(2.59) 


(2.60) 


n'  =  [f“^af  -  f"^(af  -  fj)  j ? 

(2.61) 

or. 

■?'  =  [f"^af  -  f“^af  +  f“^fj]  V 

(2.62) 

and  finally. 

"v*  =  J"? 

(2.63) 

This  represents  a  much  easier  system  to  deal  with  since  J  is 
a  constant  matrix.  If  all  the  characteristic  exponents  are 
real  and  distinct,  J  is  set  to  have  elements  only  on  the  main 
diagonal.  If  this  is  the  case,  the  system  described  by  equa¬ 
tion  (2.63)  reduces  to  n  uncoupled  ordinary  differential 
equations,  hence  the  new  n  variables  are  called  the  modal 
variables.  Each  represents  a  different  mode.  The  stabil- 
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ity  of  one  mode  is  independent  of  the  other  modes.  In  the 


state  variables  each  mode  n ^  participates  in  the  overall 
motion  of  the  system.  If  the  characteristic  exponents  are 
complex,  J  is  put  in  block  diagonal  form,  with  the  imaginary 
parts  on  the  off  diagonal.  In  this  case  two  of  the  modes 
become  coupled  and  the  differential  equations  become  more 
involved. 

The  feedback  control  term  added  to  the  basic  system, 
changes  the  mathmatical  model  to  look  like; 


x'  =  Ax  +  Bu 


(2.64) 


where  B  is  a  feedback  application  matrix  and  u  is  the  control 
vector.  Expressing  this  equation  in  modal  variables  yields; 


n*  =  Jn  +  f"^Bu  (2.65) 

The  B  matrix  is  used  to  show  which  variables  the  control 
is  applied.  The  control  matrix  B  has  the  general  form; 


'  c 

0  ... 

0 

0 

0  ... 

0 

hi 

h2*  •  • 

hn 

■  ^21 

h2*  •  • 

^2n 

(2.66) 


where  the  number  of  columns  n  is  determined  by  the  size  of  the 
control  vector  u.  The  elements  b^j  are  determined  by  the 
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physical  implementation  of  the  control.  The  control  matrix 
B  has  four  rows  because  the  state  variable  has  four  elements 


or 


-1  —  ‘ 

F  Bu  = 

I 

f, 
«« 


where  f^j  represents  the  elements  of  the  F  inverse  matrix. 
Further  multiplication  and  insertion  into  equation  (2,65), 
where 


13 

+ 

^14 

23 

+ 

^24 

-  1 

33 

+ 

^34 

[Ki  ^2  "^3  ^4  1’' 

43 

+ 

^44  _ 

(2.70b) 

n '  = 


^13  ^i4 

,  gives; 

+X^ 

K2R1 

K3R1 

K^Rl 

K1R2 

^2^2  ^2 

^3''2 

K4R2 

^2^3 

^4^3 

K1R4 

K2R4 

^3^4 

^4^4 

V 

(2.71) 


This  represents  the  controlled  system  using  scaler  control. 
Suppose  one  mode  is  unstable,  for  example  is  positive, 
By  choosing  Kg  =  =  0,  equation  (2.71)  reduces  to; 


V '  = 


^1^^13 

+ 

^14)  ♦  *'1 

0 

0 

0 

^1^^23 

+ 

^24^ 

X  2 

0 

0 

^1^^33 

+ 

^34) 

0 

X 

3 

0 

‘'l^^43 

+ 

^44) 

0 

0 

^4 

(2.72) 


Thus  the  three  stable  modes  remain  the  same,  and  the  unstable 
mode  can  be  changed  by  proper  selection  of  K^(r), 

The  differential  equation  for  the  first  mode  is; 

Vl  =  1  Ki(0(fi3(r)  +  fi4(0)  +  (2.73) 


•rnimiimmmmbi 
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Using  an  integrating  factor, 

0 

the  solution  of  (2.73)  is  given  by; 


T 

M  =  exp  \-J  +  K^(t)(f^3(t)  +  f^4(t)))dt|  (2.74) 


’’l  ~  ®^P|/  ^  ^  R('t))d't| 


(2.75) 


where  here  ri  is  the  initial  value  of  »>.»  C  is  the  constant 
part  of  the  exponential  of  equation  (2.74),  and  the  term 
R(t)  is  the  time  varying  part.  Now  by  letting  be  a  constant, 
and  expressing  the  components  of  F  inverse  in  Fourier  series 
form,  with  a  constant  and  a  time  periodic  part,  equation  (2.75) 
reduces  to; 

1,  =  (2.76) 

^  -^0 

with 


C  =  +aK^ 


(2.77) 


where  a  is  the  constant  part  of  the  Fourier  series.  Hence 
by  appropriately  selected  K  such  that  C  is  less  than  zero 
the  stability  of  mode  one  is  guaranteed. 

Now  suppose  there  is  no  constant  part  of  the  Fourier 
series,  or  that  it  is  too  small  to  be  realistically  fed 
back  by  a  constant  K.  In  this  case  a  time  varying  gain  can 
be  used  to  control  the  system.  The  integrating  factor  of 
equation  (2.74)  shows  a  K(t)  times  a  Fourier  series.  Suppose 
the  first  cosine  tern  of  the  Fourier  series  is  large  enough 
to  modify  the  mode.  A  constant  part  can  be  created  by  making 
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the  gain; 


K  (t)  =  k,  cos 
J.  -Lq  ^ 


(2.78) 


vrtiere  K.  is  constant.  V/hen  this  is  multiplied  by  the 
■^0 

Fourier  series  representation  of  (f^^g  +  ^^4)  gives; 


Ki(r)(f^3  + 


4tt 

flfgCOS  -Y“  +  . 


,,  2tt 

K.  cos 

■^0 


f  2ir  T 

l“o  — 


.  .  +  iSj^sin  +  jS^sin  +  ...  j 

(2.79) 


or 

Ki(r)(f^3  .  .  a^K.^oos  .  .^K,^oos^(S^)  * 

cos  cos  +  ...  (2.80) 

rearranging  the  second  term  in  equation  (2.80); 


a^K^^cos^(^)  =  1  +  cos  ^) 

(2.81) 

Therefore  the  constant  part  of  equation  (2.75)  becomes; 


Since  the  rest  of  the  time  periodic  terms  integrated  over 

the  period  still  equals  zero,  the  stability  can  again  be 

insured  by  selecting  the  appropriate  K  to  be  used  with  the 

0 

time  periodic  gain.  Of  course  any  sine  or  cosine  function 


25 


can  be  used  depending  on  the  Fourier  coefficient  the  designer 
wishes  to  use.  The  form  of  equation  (2.82)  remains  the  same. 

Now  suppose  there  are  two  unstable  modes  in  equation  (2.71). 
A  scaler  control  may  again  be  used  but  the  problem  becomes 
more  involved  involved.  Suppose  further  that  the  unstable 
characteristic  exponents  are  complex  conjugates.  Applying 
the  scaler  control  equation  (2.71)  becomes; 


+ 

^14) 

^2^^13 

+ 

^14)  +  « 

0 

0 

+ 

^24^  -  w 

*^2^^23 

+ 

0 

0 

7 

+ 

^34) 

^2^^33 

+ 

^34^ 

0 

/i<^43 

+ 

^44^ 

^2 ^^43 

4* 

^44) 

0 

(2.83) 


where  9  is  the  real  part  of  the  complex  conjugate  and  w  is 
the  imaginary  part.  Again  the  stability  of  modes  and 
are  unaffected.  The  and  ri^  differential  equations  are; 

1'  .  (Ki(fi3  *  ♦  .)  -i  *  (KjCfia  +  +»)  ->2 

(2.84) 

^2  ~  ^^1^^23  ^24^  ”  *^1  ^^2^^23  ^24^  +  '>  n  2 

Uhfortxinately  this  set  of  equations  does  not  lend  itself  to 
the  same  method  of  analysis  as  the  last  case.  That  is,  an 
integrating  factor  has  not  been  found  to  solve  this  set  of 
equations  and  thus  determine  what  values  of  and  Kg  will 
guarantee  stability.  However  equation  (2.84)  can  be  written 
as  a  periodic  linear  system; 

=  X(t)  7  (2.85) 
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and  then  a  theorem  from  Floquet  theory,  which  states  (3;236), 

=  ®*pj  jx(s)  I  ds  j  (2.86) 

can  offer  some  help.  Here  are  the  characteristic  multi¬ 
pliers  and  tr(X)  is  the  trace  of  the  X  matrix.  Taking  the 
log  of  both  sides  gives; 

T 

InCPiPg)  =  J  tr  I  X(s)  j  ds  (2.87) 

from  equation  (2.84)  the  trace  of  the  matrix  is; 

tr  |x(s)  I  =  +  fg^)  +  2«r  (2.88) 

or  breaking  it  up  into  constant  and  periodic  parts, 

tr  |x(s)|  *  +  Kgttgj  +  2*^)  +  R(’’)  (2.89) 

where  the  second  subscript  on  the  a's  represents  the  coef¬ 
ficient  of  the  Fourier  series  term  yielding  a  constant  part. 
These  are  determined  by  the  form  of  the  gains  and  K^, 
Integrating  the  periodic  part  over  one  period  reduces  to 
zero  while  the  constant  part  leaves; 

ln(pjP2)  =  +  Kgttgj  +  2<r)T  (2.90) 

or  by  rearranging, 

|>ln(Pj^P2)  *  +  Y^nPg  =  +  ^2*21  (2.91) 

where  1/T  In  P^^  is  recognized  as  the  characteristic  exponent 
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9 


hence 


Xi  +  Xg  = 


+  2<r 


(2.92) 


Equation  (2.91)  does  not  predict  the  value  of  the  indivd- 
ual  characteristic  exponents  X^  and  X^,  but  rather,  only 
their  sum.  Naturally  their  sum  needs  to  be  negative  for 
stability,  but  making  their  sum  negative  does  not  guarantee 
that  both  values  will  be  negative.  It  is  deireable  to  be 
able  to  predict  what  values  of  and  will  guarantee  both 
modes  becoming  stable,  but  as  of  yet  this  is  not  possible. 

The  best  that  can  be  done,  using  scaler  control,  is  to  use 
the  method  described, to  help  search  for  gains  that  will 
result  in  a  stable  system. 

There  is  a  method  that  will  predict  the  value  of  the 
controlled  characteristic  exponents.  This  method  uses  a 
multiple  input  controller,  which  enlarges  the  B  and  K  matrices 
of  equation  (2.70).  The  following  is  a  demonstration  of  a 
multiple  input  control  system  used  to  control  two  unstable 
modes.  First  let  the  B  matrix  be. 


and  the  K  matrix; 


K 


fo  0 


B  = 


0 

1 


0 

0 


[0 


1 


11 

Ki2 

•^13 

^14 

21 

'^22 

^23 

CM 
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(2.93) 


(2.94) 


K 


Hence , 


■  0 

0 

0 

0 

1 

i 

0 

0 

0 

0 

F 

BKn  = 

F”^ 

■ 

^11 

^12 

^13 

^14 

K22 

''23 

^24 

Therefore 

the  contolled  i 

system 

becomes , 

■^11^13 

^21^14 

Ki2 

^13  " 

^22^ 

TT 1 

^11^23 

^21^24 

^12 

^23  " 

K22f 

V '  = 

^11^33 

*  ^21^34 

^2 

^33 

^22^: 

.^11^43 

^  ^21 ^44 

Ki2 

^3  " 

^'22". 

+  Ul 

+  ff- 


(2.95) 


0  0 
0  0 


^3  ° 


-  ”4 

(2.96) 

where  setting  =  ^23  =  =  K24  =  0  gives  the  above 

result.  Since  the  stability  of  the  and  modes  is  not 
affected  by  the  control,  the  stability  of  the  controlled  X^ 
and  Xg  is  determined  by  the  solution  to; 


^2 


K 


11^13  ^21^14  ^12^13  ^22^14 

^11^23  ^21^24  ""  ^12^23  ^22^24 


(2.97) 


To  determine  the  values  of  the  K's  that  will  guarantee  stab¬ 
ility  of  the  system,  first  set  the  off  diagonal  terms  to  zero 
This  vincouples  the  two  modes  and  results  in. 


^11^23  *  ^21^24  “  " 
^12^13  ^22^14  =  '  " 


(2.96) 


solving  for  and  K^2  gives. 
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Next,  set  the  diagonal  terms  to  desired  constants  'a'  and 


'b'  , 


V  +  +  ^^21^14  ~  ^ 

a  +  Kj_2-23  +  ^22^24 


(2.100) 


solving  for  and  K22  yields; 


K 


K 


21 


22 


(a 

(t,  -  .)  -  Kj2f23]/f 


14 

24 


(2.101) 


Using  equation  (2.99)  in  equation  (2.101)  gives  the  K's  as; 


KjjCr)  . 


-  (a  -  «r  ) 


"■24 


■P  "P  P 

14  23  -  13'’24 


K,3(r)  = 


-“^24  -  (t>  -  <r 

f  f  -  f  f 
13  24  -^14^23 


(2.102) 


(a  -  <r  )f23  - 
^14^23  "  ^13^24 


K22(’)  - 


(b  -  «r)f^3  -H  wf23 
^13^24  "  ^14^23 


These  K's  are  time  varying  and  can  be  expressed  in  Fourier 
series  form.  Thus  the  K  matrix  with  the  above  values  of  K 
and  the  B  matrix  of  equation  (2.93)  gives  the  final  system; 


(2.103) 


where  a  and  b  are  selected  to  be  negative  and  thus  stable. 


"  fni'Mitrifia  111(1111 
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CHAPTER  3 


RESULTS 

In  this  chapter  the  results  of  various  examples  using 
the  theory  from  the  last  chapter  will  be  presented  and 
explained.  These  examples  were  done  using  various  computer 
programs  to  solve  for  the  equations  given  in  the  theory 
chapter.  These  programs  also  assisted  in  selecting  the  con¬ 
trol  system  to  be  used,  implementing  the  control,  and  simula¬ 
ting  the  results.  Presented  first  is  a  listing  and  classifi¬ 
cation  of  the  stability  parameters  for  various  uncontrolled 
satellite  configurations.  Next  the  results  of  using  scaler 
control  to  control  a  system  with  two  unstable  modes,  is 
presented.  Finally  the  case  of  two  unstable  modes  is  solved 
using  the  multiple  input  control. 

PARAMETER  SPACE  K  vS  « 

The  basic  system  was  expressed  in  the  linear  form; 

x'=A(»-)x  (3.1) 

where  A(r)  is  a  periodic  matrix  whose  elements  are  given  by 
equations  (2.26).  These  elements  are  functions  of  the  param¬ 
eters  K,  a,  and  e.  Hence  the  stability  is  uniquely  determined 
by  the  given  values  of  these  three  parameters.  In  the  article 
by  Kane  and  Barba  (4;405)  a  parameter  space  is  shown  of  K 
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verses  a  for  various  values  of  eccentricity.  In  this  study 
a  plot  of  K  verses  o  was  made  using  an  eccentricity  of  0.5. 

The  plot  originates  from  a  listing  of  the  characteristic 
exponents  for  each  pair  of  K  and  a  coordinates,  and  also  a 
classification  of  the  stability  represented  by  those  char¬ 
acteristic  exponents.  These  results  are  listed  in  Table  3.1. 
The  system  is  stable  if  all  real  parts  of  the  characteristic 
exponents  are  less  than  zero.  There  are  many  cases  v/here  the 
real  part  is  equal  to  zero,  or  very  close  to  zero,  these  cases 
are  considered  stable  since  this  represents  motion  that  is 
not  increasing  without  bound.  The  zero  real  part  is  due  to 
the  conservative  nature  of  the  system.  When  the  characteristic 
exponents  are  complex,  the  motion  will  be  oscillatory  with  the 
freqency  of  oscillation  dependent  on  the  size  of  the  real  and 
imaginary  parts.  The  damping  of  the  system  is  the  negative 
of  the  real  part  divided  by  the  natural  freqency  of  the  system. 
Hence  if  the  real  part  is  negative  the  damping  is  positive  and 
the  motion  dies  out  in  time.  If  the  real  part  is  positive  the 
damping  is  negative  and  the  motion  increases  with  timie.  The 
size  of  the  real  part  determines  how  fast  the  m.otion  either 
increases  or  decreases.  For  example,  a  large  positive  real 
part  causes  the  motion  to  increase  faster  than  a  small  posi¬ 
tive  real  part. 
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TABLE  3.1  -  CHARACTERISTIC  EXPONENTS  FOR  K  VS 


0.0000  +/-  JO. 1215  0.0000  +/-  JO. 1099  0.0000  +/-  JO. 0000  0.0000  +/-  JO. 0000 
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5384  0.3735  -0.5384  -0.3735  1.2424  0.6066  -0.6066  -1.2424 


Some  paterns  from  Table  3.1  can  be  recognized  as  to  the 
stability  dependence  on  the  parameters  K  and  a  .  Specifically 
the  real  part  of  the  characteristic  exponents  approaches 
zero  as  the  K  approaches  zero,  from  either  the  positive  or 
negative  direction.  This  means  the  instabilities  become  less 
severe.  When  K  equals  zero  the  satellite  has  the  same  moments 
of  Inertia  in  all  three  axes,  (i.e.  J=1  in  equation  (2.15)). 
When  K  equals  -1.0  the  stability  is  independent  of  the  spin 
rate.  This  case  corresponds  to  J  equalling  zero,  which  is 
not  realistically  possible.  When  J  is  nearly  zero  the  mass 
of  the  satellite  is  distributed  very  close  to  the  satellite 
spin  axis.  For  other  values  of  K,  the  larger  the  spin  para¬ 
meter,  either  positive  or  negative,  the  instabilities  become 
smaller. 

From  Table  3,1  many  of  the  characteristic  exponents  are 
in  similar  form.  The  only  cases  that  are  stable  are  those 
where  the  real  parts  all  equal  zero.  The  most  prevalent 
unstable  case  is  the  one  Involving  two  pair  of  complex  con¬ 
jugates;  one  set  with  a  negative  real  part  (stable),  the 
other  set  with  positive  real  part  (lanstable).  Another  type 
of  Instability  is  the  set  with  a  stable  set  of  complex  con¬ 
jugates  and  two  purely  real  characteristic  exponents;  one 
stable  and  one  unstable.  An  example  of  this  case  is  when 
K  =  0.8,  and  a  =  -1,0.  The  third  type  of  instability  is  when 
all  the  characteristic  exponents  are  real  with  two  of  then 
less  than  zero,  and  two  of  them  greater  than  zero.  A  few  of 
these  cases  have  been  difficult  to  work  with  due  to  the 
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numerical  problems  in  obtaining  a  periodic  solution  matrix. 
Therefore  the  characteristic  exponents  in  some  of  these  third 
cases  may  be  erroneous. 

These  different  cases  are  summarized  on  Figure  3.1,  where 
the  'o'  represents  a  stable  case,  'S'  means  single  unstable 
mode,  'C  means  complex  conjugates  unstable  case,  and  'D' 
represents  the  dual  real  unstable  case. 


FIGURE  3.1  -  PARAMETER  SPACE  K  VS 


CASE  STUDY  UNCONTROLLED 

Yeakel  presented  the  case  of  one  iinstable  mode  along  with 
a  technique  to  control  such  a  case.  A  scaler  control  is 
adequate  for  changing  this  one  unstable  characteristic 
exponent  to  a  negative  value  with  predictable  results.  From 
the  parameter  space  of  Figure  3.1  it  is  obvious  that  the  major¬ 
ity  of  unstable  cases  involve  two  unstable  modes.  As  explained 
in  the  theory  chapter  a  scaler  control  can  not  predict  v/here 
these  two  unstable  characteristic  exponents  will  move. 

The  stability  cannot  be  insured  using  a  scaler  control. 

The  only  method  to  date  has  been  to  use  the  relationship  be¬ 
tween  the  sum  of  the  characteristic  exponents  and  the  integral 
of  the  trace  of  A(’')  in  searching  for  an  acceptable  set  of 
gains  to  be  fed  back.  The  case  study  being  looked  at  has 
parameters  K  =  0.7  and  a  =  1.0,  (e  =  0.5).  This  case  represents 
the  most  common  type  of  instability,  but  the  methods  used  to 
control  this  case  along  with  the  results,  are  applicable  to 
other  similar  cases,  and  to  cases  where  the  unstable  character¬ 
istic  exponents  are  purely  real.  The  characteristic  exponents 
for  this  case  are,  (uncontrolled); 

X  =  0.0177  +/-  JO. 3667  unstable 
X„  .  =  -0.0177  +/-  JO. 3667  stable 

The  motion  of  the  uncontrolled  system  is  shown  in  Figures 
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(3.2)  thru  (3.10).  Figures  (3.2)  thru  (3.5)  show  the  motion 
of  the  state  variables,  defined  in  equation  (2.20).  Figure 
(3.6)  shows  the  motion  of  the  angle  <j>  ,  which  is  the  angle 
between  the  satellite  spin  axis  and  the  orbit  normal,  and  is 
defined  by; 


^  =  arccos(  cos6j^  00362) 


(3.2) 


Figures  (3.7)  thru  (3.10) 
ables  defined  by  equation 
can  be  predicted  directly 
The  state  variable  motion 


shows  the  motion  of  the  modal  vari- 
(2.57).  The  modf'*  variable  motion 
from  the  characteristic  exponents, 
is  a  combination  of  all  the  modes. 


I 
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FIGURE  3,2  -  Uncontrolled  x,  Response 


FIGURE  3.8  -  Uncontrolled  Second  Mode  Response 
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CASE  STUDY  SCALER  CONTROL 


From  Floquet  theory  of  chapter  2,  a  solution  to  the 
linear  system  of  equation  (3.1)  is  given  by; 

^  (r)  =  F(r)e'^"  (3.3) 


where  J  is  given  by, 


0.0177 

0.3667 

0.0000 

0.0000 

-0.3667 

0.0177 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0177 

0.3667 

0.0000 

0.0000 

-0.3667 

-0.0177 

The  periodic  F  matrix  was  found  by  computer  programs  which 
expressed  it  in  Fourier  series  form.  Each  element  of  this 
4x4  matrix  has  31  cosine  terms  and  29  sine  terms  representing 
its  function.  A  typical  element  is, 


flj^(r)  =  0.00  -  0.3239  cos(2irT/T)  +  0.00  +  0.0688  cos(6irr/T) 
+  ...  +  0.0004  cos(60irT/T)  -  0.1259  sin(2irT/T) 

+  0.00  -  0.0326  sln(6irr/T)  +  ...  + 

0.00001  sin(58»’r/T) 

where  T  is  the  period  of  the  orbit. 

Using  scaler  control  puts  equation  (2.65)  into  the  form  of 
equation  (2.83).  Controlling  the  unstable  modes  reduces 
equation  (2.83)  to  the  2x2  system  of  equation  (2.84).  The 
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trace  of  the  matrix  of  this  2x2  system,  used  in  equation 
(2.86),  is  given  by; 


f  ) 
14^ 


+  K2(f23 


^24^ 


(3.5) 


where  the  f  values  in  equation  (3.5)  are  the  Fourier  coeffic¬ 
ients  of  the  terms  which  yield  constant  parts.  For  this  case 
the  zeroth  cosine  terms  (the  constant  parts)  of  the  Fourier 
series  are  too  small  to  be  adequately  amplified  for  control. 
However  the  first  cosine  coefficients  are; 


fl3  =  0.0  -  1.0027  cos(2irT/T) 
f .  .  =  0.0  +  0.2689  cos(2rT/T) 

(3.6) 

^23  “  0.2919  cos(2jrr/T) 

fg^  =  0.0  +  0.1128  cos(2Tr/T) 

Making  the  gains  time  periodic  functions  such  that. 


K,  (’■)  =  G- cos(  2irr/T) 

^  ^  (3.7) 

K2(’-)  =  G2COs(2irr/T) 

vrtiere  the  G's  are  the  constant  gain  values  yet  to  be  deter¬ 
mined.  Now  using  equation  (3.7)  and  (3.6)  in  equation  (3.5) 
the  constant  terms  become; 


Q  =  +  Xg  =  2(0.0177)  +  Gj^(-0.3669)  +  02(0. 2024) 

(3.8) 
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Rearranging  the  terms  in  equation  (3.8)  to  express  the  gains 
in  a  ratio  gives; 


n  =  0.0354  +  G^(-0.3669  +  0.2024(02/0^))  (3.9) 

The  form  of  equation  (3.9)  allows  a  specific  i^ain  ratio  to  be 
selected  and  then  a  value  of  0^  to  be  varied  to  achieve  desired 
values  of  12,  the  sum  of  the  controlled  characteristic  expo¬ 
nents.  As  an  examnle.  for  the  ratio  of  R  =  -3.0,  if  0^  =1.0 
then  G2  =  -3,0.  From  equation  (3.9)  12  becomes, 

fl  =  0.0354  +  (1 .0) (-0.3669  +  0 . 2024( -3 . 0) )  =  -0.9387 

Using  Floquet  theory  the  new  characteristic  exponents  of  the 
controlled  system  are, 


=  -0.6883  X2  =  -0.2485 

where  the  sum  of  these  characteristic  exponents  is  -0.9368. 
These  values  show  the  system  has  been  stabilized,  and  since 
they  are  purely  real  there  is  no  oscillation  in  these  modes. 
There  is  still  oscillation  in  the  other  set  of  stable  char¬ 
acteristic  exponents,  hence  the  angle  ^  oscillates  as  shown 
in  Figure  (3.11).  For  the  same  ratio  R  =  -3.0,  if  0^  =  0.6 
then  Gp  =  -1.8  and  the  predicted  sum  would  be  -0.5491  and 
the  new  characteristic  exponents  are, 

X  =  -0.2742  +/-  JO. 1432 

1  t  2 
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FIGURE 


where  the  sum  is  -0.5484.  The  plot  of  ^  for  this  combination 
of  gains  is  shown  in  Figure  (3.12).  Here  there  are  more 
oscillations  due  to  the  imaginary  part  of  the  new  character¬ 
istic  exponents.  These  two  cases  both  exhibit  stable  behavior. 
However,  in  general  for  cases  such  as  these,  stability  is  not 
assured  and  each  case  must  be  evaluated. 

By  varying  the  for  specific  ratios,  a  locus  of  character¬ 
istic  exponents  can  be  plotted.  For  different  ratios  these 
plots  represent  a  family  of  curves.  These  are  shown  for  the 
case  being  studied  in  Figures  (3.13)  thru  (3.18).  Note  that 
all  curves  originate  at  the  unstable  uncontrolled  character¬ 
istic  exponents.  As  the  gain  is  increased  from  zero,  the 
locus  cross  the  imaginary  axis  and  become  stable.  This  cross¬ 
over  occurs  when  the  trace  of  the  2x2  system  is  equal  to  zero. 
The  locus  then  breaks  in  at  some  point  on  the  real  axis  and 
splits,  one  branch  heads  toward  negative  infinity,  and  the 
other  branch  heads  toward  the  origin.  This  break  in  point  is 
very  important  because  it  can  be  thought  of  as  the  most  stable 
placement  for  the  controlled  characteristic  exponents,  as  trill 
be  seen.  From  Figure  (3.13)  the  break  in  point  for  a  ratio 
of  0.5  occurs  around  -0.02,  and  from  Figure  (3.16),  the  ratio 
of  -10.0,  the  break  in  point  is  around  -0.35.  If  the  charact¬ 
eristic  exponents  lie  anywhere  in  the  left  half  plane,  they 
represent  a  stable  configuration.  Obviously  the  closer  they 
are  to  the  real  axis,  the  less  the  oscillation,  and  the  further 
they  are  from  the  imaginary  axis  the  quicker  the  motion  dies 
out.  The  break  in  point  is  the  point  that  is  furthest  from  the 
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FIGURE  3,13  -  Locus  of  Characteristic 
Exponents  for  R  =  0.5 


e  :  Gj^  =  0.0 

A  :  Gj^=  0,25 
*:  Gj^*  0,5 
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FIGURE  3.14  -  Locus  of  Characteristic 
Exponents  for  R  =  -0.5 


©:  G^=  0.0 

a:  G^*  0.4 

*  ;  Gj^=  0,6 


-0.1 


FIGURE  3.15  -  Locus  of  Characteristic 


Exponents  for  R  =  -3.0 


©;  Gj^=  0.0 

A:  Gj^=  0.5 
»:  Gj«  0.77 


-0.4 


-0.3 


-0.2 


-0.1 
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imaginary  axis  and  closest  to  the  real  axis  for  both  sets  of 
controlled  characteristic  exponents. 

CASE  STUDY  MULTIPLE  INPUT  CONTROL 

The  scaler  control  used  to  stabilize  a  system  with  two 
unstable  modes  is  desireable  due  to  the  simplicity  of  the  feed¬ 
back  gains  used.  However  since  there  is  no  way  to  accurately 
predict  the  placement  of  the  controlled  characteristic  expon¬ 
ents,  another  method  is  needed.  The  multiple  input  control 
system  provides  more  predictable  results.  Again  controlling 
the  unstable  modes  results  in  the  system  described  by  equation 
(2.97),  The  time  varying  gains  were  found  using  equations 
(2.102).  These  gains  are  expressed  as  Fourier  series.  The 
new  characteristic  exponents  were  chosen  to  be; 

=  -0.75  Xg  =  -0.50  (3.10) 

A  computer  program  found  the  Fourier  series  for  the  gains 
Kin  Kig,  ^21’  ^22  would  result  in  the  characteristic 

exponents  of  equation  (3.10).  An  example  of  the  time  varying' 
gains  is; 

Kii(»')  =  -1.243  cos(2tt/T)  +  0,625  cos(6»-r/T)  +  ... 

-2.742  sin(2*'r/T)  -  0.111  sin(6jrr/T)  (3.11) 

where  T  is  the  period  of  the  orbit.  A  summary  of  the  gain 
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Fourier  coefficients  for  this  case  is  listed  in  Table  (3.2). 
The  ga.ins  given  in  this  case  yield  the  nev/  characteristic 
exponents  of; 

=  -0.  7685  X^  =  -0.4777 

where  the  error  from  the  predicted  values,  is  due  to  the  num- 
ericel  inaccuracies  of  the  computer.  More  accurate  results 
v/oulcl  have  been  achieved  if  the  Fourier  series  had  been  taken 
to  more  coefficients.  The  <t>  plot  for  this  set  of  gains  is 
given  in  Figure  (3.19). 

I’ov;  if  the  new  characteristic  exponents  were  desired  to  be 

at ; 

X  ^  =  -2.00  =  -1 .50 

the  Fourier  series  of  the  gains  would  be  as  listed  in  Table 
(3.3),  and  the  resulting  characteristic  exponents  would  be  at 

X^  =  -2.235  Xg  =  -1.319 

The  0  plot  for  these  gains  is  shown  in  Figure  (3.20). 


table  3.2  -  GAIN  FUNCTION  FOR  MULTIPLE  CONTROL  CASE  1 


10  10 

K.  .(t)  =  2  A„  cos(2(n-l)Tr/T)  +  2  sin(2mTr/T) 

n=l  in=i 


^11 

^12 

^21 

K22 

cosine 
terms ; 

^1 

0.8419 

E-7 

0.6213 

E-7 

0.1670 

E-7 

0.3438 

E-8 

Ag 

-1.2429 

-0.6862 

-0.5508 

E-1 

-1.1693 

^3 

-0.5849 

E-7 

-0.1414 

E-7 

0.6755 

E-7 

0.1566 

E-7 

^4 

0.6248 

1.3368 

1.4945 

-0.2849 

^5 

-0.5868 

E-7 

-0.4116 

E-7 

-0.8632 

E-7 

-0.9058 

E-7 

^6 

1.4044 

1.6458 

1.6669 

-0.2662 

A? 

-0.1 008 

E-6 

-0.5761 

E-7 

0.3726 

E-7 

0.5176 

E-8 

^8 

-1.2509 

-0.3723 

0.5329 

-0.4886 

^9 

0.1498 

E-6 

0.1258 

E-6 

0.2521 

E-7 

0.2693 

E-8 

^10 

-2.1042 

-0.6031 

-0.3944 

-0.6961 

sine 

terms; 

®1 

-2.7417 

-1.6274 

-0.3833 

E-2 

0.2135 

®2 

0.1118 

E-6 

0.8622 

E-7 

0.6379 

E-7 

0.4292 

E-7 

^3 

-0.1112 

-0.1777 

-1.5967 

-1.2184 

S4 

0.2076 

E-7 

0.1988 

E-7 

0.1357 

S-7 

0.1058 

E-7 

-0.3721 

0.2358 

2.0491 

0.9833 

®6 

-0.1006 

E-6 

-0.9000 

E-7 

-0.5690 

E-7 

-0.4547 

E-7 

®7 

3.2091 

1.5700 

0.2321 

E-1 

-0.1490 

®8 

-0.5102 

E-7 

-0.2927 

E-7 

-0.2555 

E-7 

-0.2070 

E-7 

S9 

-2.3053 

-0.8608 

0.7180 

-0.3420 

E-1 

Bi  n 

0.9814 

E-7 

0.7360 

E-7 

0.5153 

E-7 

0.4196 

E-7 

TABLE  3.3  -  GAIN  FUNCTION  FOR  MULTIPLE  CONTROL  CASE  2 


10  10 


K.  .(t)  =  2  A„  cos(2(n-l)*-»-/T)  +  2 

ij  _  ^  n 


i  J 

n= 

=1 

n 

m=l 

cosine 
terms ; 

^11 

^12 

^21 

^1 

0.2403 

E-6 

0.1500 

E-6 

0.4410 

E-7 

^2 

-3.4555 

-1.5559 

-0.6102 

^3 

-0.1553 

E-6 

-0.2150 

E-7 

0.1791 

E-6 

^4 

2.1350 

3.6081 

3.7095 

^5 

-0.1667 

E-6 

-0.9843 

E-7 

-0.2572 

E-6 

^6 

4.2537 

4.2380 

4.1580 

A? 

-0.2811 

E-6 

-0.1318 

E-6 

0.9740 

E-7 

^8 

-3.3500 

-0.6592 

1.1674 

^9 

0.4338 

E-6 

0.3102 

E-6 

0.6560 

E-7 

^10 

-5.6257 

-1.0427 

-1.3700 

sine 
terms ; 

-7.6279 

-3.4714 

0.7585 

E-1 

0.3208 

E-6 

0.2097 

E-6 

0.1804 

E-6 

S3 

-0.3559 

-0.4705 

-4.5742 

S4 

0.6108 

E-7 

0.4998 

E-7 

0.3896 

E-7 

S5 

-0.8576 

0.7912 

5.6375 

S6 

-0.2935 

E-6 

-0.2241 

E-6 

-0.1689 

E-6 

S7 

8.8409 

3.4427 

-0.3819 

E-3 

Se 

-0.1422 

E-6 

-0.6702 

E-7 

-0.7368 

E-7 

S9 

-6.2437 

-1.7123 

1.8233 

0.2806 

E-6 

0.1782 

E-6 

0.1487 

E-6 

1  sin(  2niirr/T) 
m 


0.4434  E-8 
-3.3138 
0.2293  E-7 
-1.2916 
-0.2303  E-6 
-1.2936 
0.2793  E-8 
-1.5632 
-0.4223  E-9 
-2.4263 

0.6094 
0.1016  E-6 
-2.9586 
0.2578  E-7 
2.1437 
-0.1106  E-6 
-0.4319 
-0.5076  E-7 
-0.3279 
0.1030  E-6 


figure  3.20  -  Controlled  ^  Response  for  Multiple  Input  Control  Case 


CHAPTER  4 
CONCLUSIONS 


There  are  some  important  conclusions  which  can  be  drawn 
from  the  results  presented  in  the  previous  chapter.  Foremost 
is  the  contrast  between  the  scaler  control  and  the  multiple 
input  control.  The  use  of  these  two  control  schemes  in 
simultaneously  moving  two  unstable  characteristic  exponents 
was  demonstrated  in  the  results  section.  The  problem  of 
moving  two  roots  simultaneously  is  of  particular  importance 
to  satellite  designers,  due  to  the  number  of  cases  where  two, 
or  more,  unstable  modes  exist.  In  this  study  the  case  where 
there  are  two  unstable  modes  was  shown  to  outnumber  the  cases 
of  one  unstable  mode.  It  would  be  overly  restrictive  to 
limit  satellite  designs  to  only  those  cases  of  zero  or  one 
unstable  mode.  In  cases  with  multiple  unstable  modes  a 
scaler  control  may  be  applied  repetitiously  to  control  a 
single  mode  at  a  time.  This  involves  applying  a  scaler  con¬ 
trol  to  the  original  system  to  take  care  of  one  mode,  then 
forming  a  new  "semicontrolled"  system  to  which  another  scaler 
control  is  applied.  This  process  is  continued  lontil  all  modes 
are  controlled  as  desired.  The  problem  here  is  that  the 
numerical  difficulties  quickly  make  this  scheme  unworkable. 

The  primary  emphasis  of  this  study  has  been  to  develop  a 
workable  scheme  to  control  a  periodic  system  with  two  unstable 
modes.  From  the  results  presented  it  is  obvious  the  two  types 
of  controllers  offer  significant  differences  in  achieving 
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this  control.  The  scaler  control  works  with  simple  workable 
feedback  gains,  but  does  not  offer  predictable  results.  The 
multiple  input  control  system  offers  predictable  results  as 
to  the  placement  of  the  new  stable  characteristic  exponents, 
but  requires  complicated  gains  which  are  time  dependent.  The 
greater  the  accuracy  desired  in  the  placement  of  the  charact¬ 
eristic  exponents,  the  more  complicated  the  gain  time  functions 
become.  These  complicated  gains  are  unavoidable  in  the  multiple 
input  controller.  Therefore  the  problem  becomes,  how  to  predict 
the  placement  of  the  characteristic  exponents  using  the  scaler 
control . 

The  only  clue  in  predicting  the  results  of  the  scaler 
controller  comes  from  the  relation  between  the  characteristic 
exponent's  sum  and  the  trace  of  the  matrix  defining  the  sys¬ 
tem.  Through  the  use  of  the  theorem  given  by  equation  (2.86) 
a  locus  of  controlled  characteristic  exponents  was  plotted. 

From  this  plot  there  can  be  found  a  breakin  point,  unique 
for  a  given  gain  ratio.  If  this  breakin  point  can  be  predicted 
the  problem  will  be  effectively  solved.  Also  from  these  plots 
it  is  obvious  that  some  optimum  ratio  exists  that  gives  the 
most  negative  breakin  point  for  the  system.  As  was  discussed 
the  breakin  point  can  be  thought  of  as  the  most  stable  con¬ 
trolled  point.  More  negative  values  of  one  characteristic 
exponent  can  be  achieved,  but  it  requires  the  second  charact¬ 
eristic  exponent  to  be  less  negative.  Hence  the  breakin  point 
for  the  optimum  ratio  represents  the  most  desireable  location 
of  the  controlled  characteristic  exponents.  As  of  yet  there 
is  no  clear  way  to  predict  breakin  points  or  even  the  optimum 
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ratio . 


If  a  less  specific  solution  is  desired,  that  is  one  that 
is  simply  stable  and  not  necessarily  the  most  stable,  then 
the  imaginary  axis  crossings  must  be  predicted.  For  cases 
similar  to  the  one  presented,  (  that  is  one  with  negative 
breakin  points),  the  first  imaginary  crossing,  as  the  gain 
is  increased  from  zero,  occurs  when  the  sum  computed  from 
the  trace  theorem,  is  equal  to  zero.  The  second  crossover 
occurs  as  one  branch  heads  toward  the  origin  after  breaking 
into  the  real  axis.  This  crossover  was  difficult  to  find 
using  computer  numerical  integration.  Whether  a  crossover 
actually  occurs  at  the  origin  is  still  unclear.  These  two 
crossover  points  define  the  region  where  the  system  is 
guaranteed  stable. 
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